function R = cosine_pump(obj, amp, freq, offset)
%COSINE_PUMP Summary of this function goes here
%   Detailed explanation goes here
    R = chebfun( offset + amp*cos(2*pi*freq* obj.t) );
end

